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SUMMARY 

Assuming  a  constant  circulation,  rigid,  trailing 
wake  and  using  a  lift  deficiency  function  to  represent  the 
unsteady  aerodynamic  effects,  various  methods  of  calculating 
the  airloads  on  a  helicopter  rotor  in  steady,  forward  flight 
were  developed  for  the  purpose  of  achieving:  (a)  faster 
solution  times  for  a  given  level  of  accuracy  and  (b)  a  more 
accurate  representation  of  the  lifting  surface/vortex  inter¬ 
action  for  the  case  where  the  blade  passes  close  to  the  vortex 
line  generated  by  a  previous  blade.  Since  most  of  the  solution 

,  v* 

time  is  required  to  calculate  the  induced  velocities  due  to 
the  trailing  wake,  various  approximate  methods  of  calculating 
the  induced  velocities  due  to  a  rigid,  skewed  helix  were 
developed.  Within  the  limits  of  accuracy  available  from  the 
rigid  wake  model  either  the  solution  of  Reference  (2) ,  which 
represents  the  skewed  helix  with  a  series  of  infinite  straight 
lines  placed  tangent  to  the  helix  whenever  the  helix  passes 
under  a  blade,  or  the  finite  straight  line  solution,  which 
represents  the  helix  with  a  series  of  straight  line  segments 
subtending  equal  A<^  (change  in  azimuth  angle)  ,  using  Ajzi  =  30°- 
40°  were  found  to  give  the  fastest  results.  If  the  model  is 
refined  to  include,  for  example,  a  nonrig id  wake,  then  the 
finite  straight  line  solution  with  a  of  approximately 
15°  should  be  used.  The  use  of  a  conventional  two  dimensional 
lifting  surface  theory  was  tried  for  the  lifting  surface/ 
vortex  interaction  case  and  little  improvement  resulted.  It 
is  recommended  that  the  use  of  more  sophisticated  lifting 
surface  theories  be  investigated  for  this  case. 
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SYMBOLS 
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location  of  an  element  of  vorticity  in  the  trailing 
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culated  . 

rotor  radius. 

tip  path  plane. 

velocity  perpendicular  to  blade  in  TPP. 

P 

W*/  - —  )  nondimensional  induced  velocity  due  to  a 

4TTR 

vortex  segment. 

induced  velocity  due  to  a  vortex  segment. 

Wn . 

vectors  from  P^  to  the  ends  of  the  vortex  segment. 


i 

zero  harmonic  of  flapping. 

blade  semichord. 

a  -  b. 


nondimens ional  distance  travelled  by  the  rotor  hub 
between  t ^  and  t  ^  . 


h*/R. 

perpendicular  distance  from  to  the  vortex  segment. 
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instant  at  which  the  element  of  vorticity  at  P^  was 
trailed . 

instant  at  which  the  induced  velocity  at  P^  is 
calculated . 


Zlim  '>Vj: 


chordwise  induced  velocity  distribution. 


x  -  components  of  the  distances  from  to  the 
ends  of  the  vortex  segment 

chordwise  coordinate  (Section  III)  . 


y-components  of  distances  from  P 
of  the  vortex  segment. 


'I 


to  the  ends 


z-components  of  distances  from  P^  to  the  ends 
of  the  vortex  segment. 


circulation  of  trailing  vortex  (constant) . 

angular  velocity  of  rotor,  radians/sec. 

angles  between  lines  from  P^  to  the  vortex  segment 
and  the  vortex  segment  (see  Figure  23) . 


angles  between  the  tangent  to  the  spiral  and  the 
two  half-infinite  vortices  (see  Figure  22) . 

angle  between  the  induced  velocity  vector  and  the 
perpendicular  to  the  TPP. 

angle  between  the  blade  which  trailed  the  element 
of  vorticity  and  the  blade  on  which  induced  velocity 
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is  being  calculated,  at  any  instant. 
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=  angle  in  the  TPP  between  the  local  wind  and  the 

perpendicular  to  the  blade  at  point  P^  . 

1 

=  nondimens ional  mean  inflow. 


*\  "\  =  nondimens ional  n  ^  harmonics  of  inflow. 

/xytc  )  »'ns 


=  advance  ratio  (tip  speed  ratio) . 


=  <f>-  ^  -  6  (see  Figure  14) . 


=  See  Figure  16. 

=  azimuth  angle  of  blade  trailing  the  element  of 
vorticity  at  point  P^  at  the  instant  of  trailing 
(time  t ^  )  . 

=  change  in  azimuth  angle  between  the  ends  of  a  vortex 
segment . 


=  change  in  azimuth  angle  between  the  nearest  point  and 
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and  the  spiral  (see  Figure  16) . 
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I .  INTRODUCTION 


This  report  discusses  various  methods  for  computing 
blade  airloads  which  were  developed  for  the  purpose  of  achieving: 
(a)  faster  solution  times  and  (b)  a  more  accurate  representation 
of  the  lifting  surface/vortex  interaction  for  the  case  where  the 
blade  passes  close  to  the  vortex  line  generated  by  a  previous 
blade . 


The  original  solution  for  the  induced  velocity  in  forward 
flight  due  to  the  trailing  wake  used  numerical  integration  down 
the  spiral  wake  (Reference  1) .  This  method  required  small  inter¬ 
val  sizes  (typically  7-5°  in  azimuth)  and,  hence,  large  amounts 
of  computer  time  to  get  accurate  results.  A  solution  was  also 
developed  (Reference  2)  where  the  spiral  wake  is  replaced  by  a 
set  of  infinite  straight  line  vortices  (for  which  the  induced 
velocity  is  known)  placed  tangent  to  the  spiral  at  every  point 
where  the  spiral  passes  under  a  blade  during  the  first  turn  of 
the  spiral.  This  solution  is  20  times  as  fast  as  the  numerical 
integration  solution;  however,  it  is  less  accurate.  Attempts 
were  then  made  to  develop  approximate  solutions  which  would  be 
faster  than  the  numerical  integration  solution  and  more  accurate 
than  the  approximate  solution  outlined  above.  Two  such  attempts 
are  described  below. 

Nearest  Point 


Instead  of  replacing  the  spiral  trailing  wake  by 
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infinite  straight  lines  placed  wherever  the  spiral  passed  under 
a  blade,  it  was  decided  to  locate  the  infinite  straight  lines 
wherever  the  distance  from  the  point  at  which  the  induced  velo¬ 
city  was  being  calculated  to  the  spiral  was  a  local  minimum 
or  maximum  (nearest  points) .  As  a  further  refinement,  a  double 
infinite  straight  line  solution  was  developed  which  used  two 
half  infinite  straight  line  vortices,  both  starting  at  the 
nearest  point  and  proceeding  in  nearly  opposite  directions  to 
infinity,  intersecting  the  spiral  at  points  away  from 

the  nearest  point  on  both  sides  (see  Figure  1) .  An  approximate 
lifting  surface  theory  was  also  developed  based  on  the  single 
infinite  straight  line  geometry.  Appendix  A  discusses  these 
solutions  in  more  detail.  Since  these  solutions  proved  to  be 
both  slower  and  less  accurate  than  the  finite  straight  line 
solution  discussed  below,  they  were  discarded. 

Finite  Straight  Line  (FSL) 

This  solution  replaces  the  spiral  vortex  trailing 
wake  with  a  series  of  straight  line  segments  (see  Figure  2) . 
Originally,  it  was  intended  to  take  shorter  line  segments  in 
the  vicinity  of  nearest  points  and  longer  segments  elsewhere. 
Due  to  the  considerable  amount  of  computer  time  involved  in 
finding  all  the  nearest  points,  however,  a  solution  using  a 
constant  &<f>  (change  in  azimuth  angle)  per  line  segment  proved 
to  be  faster  for  equivalent  accuracy. 

Since  FSL  is  approximately  six  times  as  fast  as  the 


numerical  integration  method  for  equivalent  accuracy/  it  has  been 
adopted  as  the  normal  method  of  calculating  induced  velocities. 

It  is  even  used  for  distorted  wake  cases  where  the  trailing  wake 
is  not  a  perfect  spiral  (skewed  helix)  but  some  more  general 
shape  (Reference  3) .  Appendix  B  discusses  FSL  in  more  detail. 
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II.  RESULTS  AND  CONCLUSIONS 

Since  the  airloads  on  the  H-34  rotor  are  available 
from  NASA  Flight  Test  data  (Reference  4) ,  the  airloads  for  a 
test  case  (ji  =  . 2,  "X  =  .025)  have  been  calculated  using  in¬ 
duced  velocities  calculated  by  the  various  approximate  methods 
outlined  above.  This  gives  a  standard  against  which  the  various 
methods  can  be  compared.  The  airload  calculations  use  a  two 
vortex  (one  at  the  tip  and  one  at  the  midspan  point)  constant 
circulation  trailing  wake  with  the  unsteady  aerodynamic  effects 
accounted  for  by  the  appropriate  lift-deficiency  function 
(Reference  2) .  Since  the  lower  harmonics  of  airload  can  be 
calculated  reasonably  accurately  by  uniform  inflow  methods/ 
the  airloads  have  been  plotted  with  0th,  1st,  and  2nd  harmonics 
extracted  to  emphasize  the  higher  harmonics  which  uniform  in¬ 
flow  methods  cannot  calculate  accurately.  Figure  1,  4,  5,  and 
6  show  the  airloads  at  the  95%  spanpoint  (  Tj  =  .95)  as  cal¬ 
culated  using  the  various  approximate  methods  of  calculating 
induced  velocity  outlined  above  compared  with  the  NASA  Flight 
Test  data  (Reference  4) .  For  comparison  the  older  infinite 
straight  line  (Reference  2)  and  numerical  integration  (Reference 
1)  methods  are  shown  in  Figures  7  and  8(a)  respectively.  From 
these  figures,  it  can  be  seen  that  although  the  various  methods 
do  not  agree  completely  with  each  other,  there  is  little  dif¬ 
ference  as  far  as  the  ability  to  predict  the  experimental  data 
is  concerned.  The  finite  straight  line  case  (Figure  6)  is  shown 
for  an  interval  size  (  )  of  15°.  A  case  with  =  7.5° 

has  been  computed  which  gives  essentially  the  same  results,  and. 
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hence,  convergence  to  the  case  A  f-*-0  has  occurred  and 
Figure  6  results  in  the  correct  solution.  Since  Figure  6  does 
not  show  complete  agreement  with  the  test  results,  the  two 
trailing  vortex,  constant  circulation,  rigid  wake  model  used 
is  inadequate,  the  ability  of  the  finite  straight  line  ap¬ 
proximation  to  represent  the  model  being  better  than  the 

ability  of  the  model  to  represent  the  actual  case.  Therefore, 
the  model  requires  further  refinement,  for  example,  as  has  been 

done  in  Reference  3  which  discusses  initial  efforts  to  improve 
the  model  by  finding  a  better  approximation  to  the  wake  geometry 
than  the  assumption  of  a  rigid  wake. 

It  is,  therefore,  concluded  that:- 

1)  If  the  two  vortex,  constant  circulation,  rigid  wake 

model  is  used,  either  the  infinite  straight  line 
approximation  of  Reference  2  (which  is  fastest)  or 
the  finite  straight  line  approximation  with  a  large 
interval  size  (  =  30°  -  40°)  should  be  used  since 

the  other  methods  discussed  above  are  all  slower  and 
yield  no  better  agreement  with  the  Flight  Test  data. 

2)  Future  work  should  concentrate  on  developing  the  model 

to  include  refinements  such  as  a  non-rigid  wake  geometry 
and  a  varying  strength  trailing  wake.  For  these  cases 
with  a  refined  model,  the  finite  straight  line  aproxi- 
mation  with  an  interval  size  (A  <j>  )  of  approxi¬ 

mately  15°  should  be  used. 
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III.  TWO  DIMENSIONAL  LIFTING  SURFACE  THEORY 


When  a  vortex  line  passes  close  to  a  rotor  blade,  as 
occurs  on  mostly  lightly  loaded  rotors,  the  lifting  line  re¬ 
presentation  of  the  rotor  blade  becomes  of  questionable  accuracy 
due  to  the  rapid  variations  of  induced  velocities  over  the 
blade  chord. 

In  the  Appendix  of  Reference  1,  a  two-dimensional 
lifting  surface  theory  is  developed.  The  induced  velocity  is 
written  as  . 

,  &  + 

-  z  f  £,  1 

where  x  =  b  cos  0  is  the  chordwise  coordinate.  It  is  then  shown 
that 

l  -  y  k  [  i  h  ( *  «,) 1  rt  (n. + i  n,H 

The  A's  can  be  calculated  by  calculating  the  induced 
velocity  at  a  number  of  chordwise  stations  (typically  7)  and 
doing  a  chordwise  harmonic  analysis.  This  process  takes  ap¬ 
proximately  7  times  as  long  as  the  equivalent  lifting  line 
procedure . 


It  was  observed  that  lifting  surface  theory  results 


only  differed  significantly  from  lifting  line  results  in  the 
vicinity  of  sharp  peaks  in  the  airloads.  A  program  was,  there¬ 
fore,  developed  which  normally  used  lifting  line  theory.  When¬ 
ever  the  value  of  the  integrand  more  than  doubled  or  halved  in 
one  integration  interval  (  A  <p  )  the  program  shifted  to  lifting 
surface  theory  and  when  the  integrand  no  longer  doubled  or 
halved  in  one  the  program  shifted  back  to  lifting  line 

theory.  This  program  was  written  using  the  method  of  calcu¬ 
lating  induced  velocities  by  numerically  integrating  over  the 
wake  which  was  developed  in  Reference  1. 

Figure  3  shows  the  comparison  of  lifting  line  and 
lifting  surface  airloads  with  harmonics  below  the  third 
eliminated  compared  with  test  data  (Reference  4) .  This  is  for 
a  rigid  wake,  4-blade,  JL  =  .2,  "X  =  .025,  and  =.95  case. 
There  is  some  effect  due  to  the  lifting  surface  theory  and  a 
marginal  improvement  in  agreement  with  the  experimental  data. 

It  is  not  felt  that  this  particular  lifting  surface  theory  is 
worth  pursuing  further.  Instead,  work  is  proceeding  on  more 
complete  lifting  surface  theories  (Reference  5) . 
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APPENDIX  A 
NEAREST  POINT 


1.  Finding  Nearest  Points 


First  it  is  necessary  to  find  a  set  of  local  nearest 
points  for  a  rigid  wake.  Consider  a  case  where  one  wishes  to 
calculate  the  induced  velocity  at  time  t^  and  at  point 
due  to  vorticity  trailed  at  time  tj_  from  the  nearest  point  P^ 

•  t 

Point  P^  is  at  a  radial  distance  R  from  the  rotor  hub  on 
a  blade  with  azimuth  angle  j/  at  t^  ,  where  R  is  the  rotor 
radius.  Point  is  at  a  radial  distance  JL R  from  the  hub 

on  a  blade  with  azimuth  angle  $  at  t^  .  Let  5  ke  the  angle 
between  the  vj  blade  and  the  L  blade  at  any  instant  (  5  =  <f> 

when  t^  =  tL  )  .  Then  t^  -  t^  =  ^-  (  ^  +  S  -  ^  )  where  Jl 
is  the  angular  velocity  of  the  rotor  (d^/dt)  .  Therefore,  the 
distance  travelled  by  the  rotor  hub  in  time  t^  -  t^  is  dR  = 
v  (t  ^  ~  ti  )»  where  V  is  the  hub  velocity.  If  JL  =  V/Jl'R  = 

advance  ratio:  (see  Figure  9) 


a  =  y/(p+  S  -tj>  ) 


The  condition  for  P ^  to  be  a  nearest  point  (or  farthest  point) 
is  that  the  line  from  Py^  to  P^  (of  length  a)  be  perpendicular 
to  the  tangent  to  the  spiral  vortex  at  P ^  .  Define  0  as  the 

angle  between  the  tangent  to  the  spiral  vortex  at  P ^  and  a  line 
perpendicular  to  the  jl  blade  at  P^  .  From  Figure  9  it  follows 
that : 
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^  sin  ^sin <f>  -  a  sin  {  <f>  -  Q)  (A.l) 

^  cos  ^  -  Jcos  j  -  a  cos  {  <j>  -  0  )  =  d  =  n  ( ^  -  <j) ) 


To  find  ©  it  is  assumed  that  the  vortex  is  trailed  parallel  to 
the  local  velocity  vector,  which  is  composed  of  two  components: 
a  velocity  JtjiR  perpendicular  to  the  blade  due  to  the  rotation 
of  the  blade  and  a  velocity  parallel  to  the  hub  motion  V  = JX  Si  R 
due  to  the  hub  motion. 


From  Figure  10  it  can  be  seen  that: 

tan  6  =  -/■  C°S  ^ -  (. 

✓6+  J*  sin  j 

This  set  of  Eqs.  (A.l),  (A. 2),  and  (A. 3)  can  be  combined  into 

one  equation: 


Ji  {  A  sin  (£  +  jx  )  (  ^  +  "S  -  <j>  )  -  JX  cos  (j/  +  cos 

+  Jt  ^  sin  ( jA  -  ^  )  =0 


(A. 4) 


which  can  be  solved  for  <p  on  a  digital  computer  by  a  Newton- 
Raphson  iteration.  However,  there  are  many  solutions  over  the 

intervals  of  interest  (e.g.,  j>  =  0  - *-4ir)  and  finding  all 

the  solutions  in  such  an  interval  takes  a  considerable  amount  of 
time . 


'■w  ^ 


«  / 


> 
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Having  solved  Eq.  (A. 4)  for  one  can  now  solve  Eq. 
(A. 3)  for  ©  and  the  Eq.  (A.l)  for  a. 


To  calculate  the  downwash  P  the  following  equation  from 
Appendix  B  may  be  used: 


For  an  infinite  straight  line  vortex  oC  = ^  =  0.  From 

Figure  11: 

2  2  2  C 

h  =  Z  +  a  ,  P  =  -  W  cos  d 

r  ,  -  2  a 

cos  d  =  a/h,  P  =  — - - - - 

Z  +  a 

Where  the  downwash  P  is,  the  induced  velocity  component  perpendi¬ 
cular  to  the  tip  path  plane  (TPP) ,  and  z,  the  nondimensional 
displacement  of  the  vortex  below  the  TPP,  is  given  by  (for  a  rigid 
wake)  : 
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z  -  >v  +  S  -  <f>  > 

For  the  near  wake  case  the  vortex  line  is  only  half¬ 
infinite,  hence  cC=  0  but  f  0.  From  Figure  12: 

A  =  TT/2  -  ©  ,  h  =  (i-^  )  cos  e 

p  _  1  4-  sin© _ 

(JL  -  )  cos  0 

Note  that  z  =  0  here,  and  hence  P  =  W. 


b.  Lifting  Surface 

In  the  Appendix  of  Reference  1,  a  lifting  surface 
theory  is  developed.  The  induced  velocity  is  written  as: 

oo 

A  A 

W  =  AT  (x  )  =  cos  n  0 

where  x  *  b  cose  is  the  chordwise  coordinate.  It  is  then 
shown  that: 

L  =  st 'srJL  ^ it 

+  V((U|R,)]. 


12 


n  "* 

and  then: 


If  fj j  ~~  and  are  neglected  compared  to 


L  -  yi  t  !  I  TT  +  YU.] 


JUy 

Noting  that  SI  -  -rf-  then: 

<31  "t 

Jtf.  if  ML 
if  dt  ' 


Also: 


V  rr  Y  era  6  i-  Jlft  Sin  f 


Hence : 


w+hM 


cos 


9  +  JlR  s 


nf)  #J. 


Figure  13  shows  the  type  of  error  involved  in  neglecting 

M, 


d 


1  ; 


and  .  This  figure  is  the  result  of  calculating  the 

downwash  at  seven  chordwise  stations,  performing  a  chordwise 
harmonic  analysis  on  the  result,  and  using  the  resulting  coef¬ 
ficients  to  calculate  L. 


Aq  is  simply  the  average  value  of  the  downwash  as  a 
function  of  x  at  any  and  .  Thus  Aq  can  be  found  by 

integrating  the  downwash  over  the  chord  and  dividing  by  the  chord. 
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fL  = 


^/zft 

_  r  za,(x)  h 

Zlk/R)  J  .1Z  +  a, *6 0 


2(? 


From  Fiqure  14: 


ct.(X)  =  -  a  +  X  sinAj 


where: 


=  ^  -  © 


Since  the  chord  is  small  relative  to  the  rate  of  change  of  z, 
one  can  assume  that  z  is  constant  over  the  chord.  Substituting 
and  integrating: 


ML  __  it  3*  _  V?. 

if  ~  5-0  if  5-o 


*  • 
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cqso  .  Uz\4-*1-'IZ*-{k)siihJ  +  ?/f)  s'>n^  ) 

•'•  if  =  Z^anN)  M^?V  Ul+ 4.a(L)  s, nO  +(%]*■  sin3-*  J 


I 


2.f^)sinAj(^S:t4-4-A1-|Za.(^)sin'J  +  l(jf  6- inzo 

I  t-  (M  ufiY'-  \  (4-i+UZ- ft '<">  +  1(t)  *" 

(*'u«,e,“1  + 18^ s,n,)  U\U$U««  +$)W«. 


'  f 2.  («■  )  s'n 


^  cosaJ  -f 


4- <31  14)  Cos  A)) 


For  the  near  wake  it  is  assumed  that  the  vortex  originates 
at  the  same  chordwise  coordinate  x  as  the  point  where  downwash 
is  calculated.  From  Figure  15: 

o4=0  yyS  =  y  +  0  j  «nJ 

Thus  for  the  near  wake  case,  use  the  far  wake  equation  multiplied 
by  J~  sm&  and  with  a.-  —  (Jr cos  0  • 

3.  Double  Infinite  Vortex 


parts 


For  this  case,  the  downwash  (P) 

(P  and  P  ) ,  one  for  each  vortex. 
4"  — 


is  calculated  in  two 
From  Figure  16: 


J  sm(  <)>  -$</>)  —  1 ^ 

J  cos  (  i  ijlty 


i  sm<j>  ~  J  sin (<j> 4- $<j>) 

JcoS  j  -h jLb  f  —  yf  COS  fay) 


'e'^> 

z  <f>  Ji -e +  . 
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From  Figure  16  it  can  be  seen  that  h  =  a  cos  Yj 


<x:+  =  t-^-a;+  o  }  o  )j&_= 

-A.  —  COS  ($-  <0]  p  -  <k  COS  Y_[\+eos(<£-,oJ) 

*  •  r 4-  -7^  I  2.  ->  \/  )~  _ V 


'+  2  +  a1 


COS2-  Yi 


^+0iZco^Y„ 


For  far  wake  P  =  P  +  P  and  for  near  wake  P  =  P, . 

■f  —  + 
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APPENDIX  B 

THE  INDUCED  VELOCITY  OF  A  VORTEX  SEGMENT 


To  calculate  induced  velocities,  use  may  be  made  at  the 
following  relationship  for  a  straight  line  vortex  segment,  derived 
from  the  Biot  Savart  Law  (Reference  6) . 


y*= 


r 


4-xh 


COS  d.  4-  COS yS) 


where  W*  is  the  dimensional  induced  velocity  at  point  P^  , 
perpendicular  to  the  plane  determined  by  the  vortex  segment 
and  point  P_  ,  P  is  the  dimensional  circulation  of  the  vortex 
segment,  h*  is  the  perpendicular  distance  from  the  vortex  segment 


to 


point  P^and  oL  and are  the  angles  between  lines 


drawn  from 


the  ends  of  the  vortex  segment  to  point  Py  and  the  vortex 
segment  (see  Figure  17) .  Now  define  the  nondimens ional  quan¬ 
tities  : 


VI  2  W 


where  R  is  the  dimensional  rotor  radius.  Therefore: 


VI  = 


Cos  d  -f  Cos 


h 


i 
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and  time  t. 


APPENDIX  C 

FINITE  STRAIGHT  LINE  (FSL) 

Consider  calculating  the  induced  velocity  at  point  P^ 

.  Point  P.,,  is  on  the  quarter  chord  line  and  at 

i  ( 

a  spanwise  distance  from  the  rotor  hub  v|R  on  the  blade  at  az¬ 
imuth  angle  (J/  at  time  .  To  calculate  the  induced  velocity 

due  to  the  entire  wake,  consider  a  general  vortex  element  at 
point  ,  which  is  inducing  part  of  the  total  induced  velocity 

at  P^  .  Let  t|  be  the  time  at  which  the  element  of  vorticity  at  ^ 
was  trailed.  Then  at  time  t^  point  P^  was  on  the  quarter 
chord  line  and  at  a  spanwise  distance  from  the  rotor  hub  JL  R 
on  the  blade  at  azimuth  angle  <f)  .  Let  the  azimuth  angle  between 
the  blade  on  which  P^  is  located  and  the  blade  which  trailed  the 
element  of  vorticity  at  at  any  instant  be  2  •  Hence,  in 
the  time  from  t ^  until  t^  the  rotor  hub  will  move  a  nondimen- 
sional  distance  d  =  S  -  I  )  *  where  yi  =  is  the 

advance  ratio.  Similarly,  the  vertical  displacement,  perpendi¬ 
cular  to  the  tip  path  plane  (TPP) ,  of  point  P ^  relative  to 
point  P^  at  time  t^  will  be,  nondimens ionally, 
where  is  the  nondimens ional  mean  inflow  through  the  rotor. 


Define  a  right-handed  axis  system  with  origin  at  P-^ 
The  y-azis  and  j  are  parallel  to  and  in  the  opposite  direction 
from  the  line  of  rotor  hub  motion,  and  are  in  the  TPP.  The 

— Jk 

x-axis  and  i  are  perpendicular  to  the  y-axis  and  are  in  the  TPP. 
The  z-axis  and  k  are  perpendicular  to  the  TPP  and  point  downv/ard 
(Figure  24).  Let  i,  j,  and  k  be  unit  vectors. 
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For  convenient  calculation  the  spiral  wake  is  now 
broken  up  into  a  series  of  straight  line  segments  of  equal 
A  (ji  (typically  15°)  ,  where  Ay(  refers  to  the  difference  in 

between  the  ends  of  the  vortex  segment.  P^  is 
defined  to  be  at  the  end  of  the  vortex  segment  nearest  along 

the  spiral.  Define  the  coordinates  of  the  ends  of  a  typical 
vortex  segment  to  be  (x,  y,  z)  and  (x^y^z^)  (Figure  18).  From 
Figure  19  and  recalling  the  values  of  z  and  d: 

X  =  Jsjn  <f>  -  rj  s'm  f  l  X,  =■  i  S\n  (^  ~  ^  ~  7  s’nf' 

)/}  ^  J  cos  (</>-&<!>)  -Jcos  if  -lyutylS- 


azimuth  angle  <f> 


where  aQ  is  the  coning  angle. 


From  Appendix  B,  the  nondimensional  induced  velocity  is: 


Cos  oC 


4.  CoS 


Define  a  and  b  as  vectors  from  P^  to  the  ends  of  the  vortex  segment 
and  c  =  a  -  b,  then  from  Figure  18: 
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(-a)-(-c) 


COS  oL 


<rt) 


♦  C 


=  \t\\c\  )  C0S/  "  ltllcl 

\  K-c)x(-t)l 

h  ®  | <x\  sinoc-  -  I I 


Define  n  as  the  unit  vector  parallel  to  the  induced  velocity  at 


,  then  from  Figure  18: 


a. 


*  l 


Y\  = 


la  xtl 


The  induced  velocity  vector  is: 


;=  V/  Y?  ~ 


X-t).(-t)  ,  (-£)< 


m 


+ 


it  i 


K-c)X(-S)l 


ftxl 

\<x  x  1 1 


\ 


From  Figure  18: 


a.=  xt+'/^+2  k  X,  t  +  +  t.k 


=  (x-x,)  t  +  (y-y,)  j  +■ 


Substituting  and  reducing  results  in: 
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Vt-  xx, - 

i  *?*  Y,V  V  Aj 

tZ 

(*,2-  X*,)l+| 

in  -  mi 

*  i 

G/VY.0 1 1  +  (xy,-x,y))<  _ 

which  can  easily  be  programmed  for  rapid  computation  on  a  digital 
computer. 


Now  the  W  for  the  various  segments  of  wake  must  be  summed 
over  <£>  between  some  limits.  The  upper  limit  is  the  upper  end  of 
the  spiral  vortex  at  time  t ^  .  After  several  turns  of  the  spiral 

P becomes  so  far  away  from  that  it  induces  a  negligible 

velocity  at  P^  and  can  be  neglected.  A  typical  lower  limit  is 
two  turns  of  the  spiral  (720°  in  <f>)  below  the  upper  limit.  After 
summing  over  <j)  summing-  over  Ji  (for  two  trailing  vortices 
£  =  1.0,  .5,  typically)  and  S  (for  4  blades  3=  0,T/2/TT"/  and 
3  T  /2)  gives  the  nondimens ional  induced  velocity  at  P., 
rigid  trailing  wake  of  constant  circulation. 


due  to  a 
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FIG.  14  LIFTING  SURFACE  GEOMETRY 
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FIG.  15  LIFTING  SURFACE  ,  NEAR  WAKE  GEOMETRY 
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FIG.  19  SPIRAL  VORTEX  -  VORTEX  SEGMENT 
GEOMETRY  PROJECTED  INTO  TPP 


